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Abstract 

The effect of free carriers on dispersion and damping of coupled phonon-plasmon modes 
is considered in the long-wave approximation. The electron and phonon scattering rate as 
well as Landau damping are taken into account. 
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I. INTRODUCTION 



After the pioneering paper by Migdal JL], the effect of electron-phonon interac- 
tions on the phonon dispersion has been permanently debated. The problem was 
how to explain that the many-body approach based on the Frohlich Hamiltonian 
gives a strong phonon renormalization. For instance, the sound velocity is renor- 
malized on the order of phonon-electron coupling A which has the value of unity. 
It follows that the phonon-electron system can be unstable. The phonon softening 
caused by the electron-phonon interaction have been discussed in many papers. All 
these results contradict to the Born-Oppenheimer conception ||] according to which 
the light electrons must follow the slow lattice vibrations. Therefore, the phonon 
renormalization should involve the small parameter of nonadiabacity \Jm/M, where 
m and M are the electron and ion masses correspondingly. 

The discrepancy was resolved by Brovman and Kagan ||. They demonstrated 
the shortcomings of the Frohlich model. Employing the adiabatic approximation, 
they found that there are two terms in the second order perturbation theory, which 
compensate each other making a result small by the nonadiabatic parameter. 

But recently, Alexandrov and Schrieffer again obtained the strong phonon 
renormalization and predicted an extremely large dispersion of optical phonons (on 
order of the Fermi velocity) because of coupling to electrons. The large phonon 
dispersion is a typical result of the theory || used the Frolich Hamiltonian. Reizer 
H stressed the importance of screening effects accompanying the longitudinal optical 
modes, but the nonphysical renormalization remained in his results. For the first 
time, the screening of the Coulomb field in optical vibrations was treated in the 
paper 0. 

Only the phonon frequencies were calculated in the works [[|, |J, 0, with no 
results available for the attenuation of optical phonons. The electron collisions with 
each other, with defects and phonons were ignored. Beside the electron collision 
rate 7, the natural phonon width r nat was also disregarded. The natural phonon 
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width is caused by the anharmonicity processes of the phonon decay into two (or 
more) phonons. Notice, that the collision processes determine the conductivity and 
the dielectric permittivity, i. e. electrodynamics of the electron-phonon system. In 
the optic range, the collision rates 7, r nat ~ \j m l Muo are small compared with a 
typical phonon frequency ujq and they give the widths of the plasmon and phonon 
resonances. Experimental studies of this resonances provide informations about iso- 
topic compositions and quality of semiconductor materials. 

Using the Boltzmann equation for electrons and the equation of motion for 
phonons we calculated || the frequency shift and the width of optical phonons in 
metals. We take into account all the mentioned features: the electron-phonon inter- 
action of different forms, the Coulomb screening, and the collision rates of electrons 
and phonons. We find that dispersion and damping of the longitudinal modes can be 
correctly described if we neglect the direct electron-phonon interaction A and keep 
only the screening and the collision rates 7, r nat . Emphasize, that in the semiclas- 
sical approximation, when the phonon momentum transferred to electrons is small 
compared with the electron momentum, the method of the Boltzmann equation is 
completely equivalent to the diagram technique. The corresponding equations may 
be formulated as equations for the electron and phonon self-energies. But in any 
case, we must properly incorporate adiabatic approximation and screening. 

Experimentally, it is convenient to investigate the free-carrier effect on phonon 
modes, varying the carrier concentration, i. e. for doped semiconductors or super- 
conductors (see, for example, the recent work [fL0|l , carried out on the HTcS com- 
pound Ndi. 86 Ceo.i4Cu04 +( 5 using IXS). Therefore, in the present paper, we extend 
the results of the paper |J to the case of small free-carrier concentrations when the 
electron plasma frequency u ep can be of the order of the phonon frequency uo, and 
the coupled phonon-plasmon modes can exist. We consider the screening as well as 
the collision effects 7 and r nat , focusing on the width of the phonon-plasmon modes. 
We simplify essentially the problem ignoring the direct electron-phonon interaction 
(terms with A) and assuming that, first, the electron system is degenerate and, sec- 
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ond, the momentum transfer k is small compared with the Fermi momentum pp. 



II. ASYMPTOTIC EXPRESSIONS FOR THE DIELECTRIC FUNCTION 

Let us find the limiting expressions of the dielectric function 

^ 2 -^ o + ,r nat Ane 2 u (v z /A p (k)) 



e(k, lo) = e 



DC 



J o; 2 -a;| +ir nat fc(l - i(j/A p {k)}) ' v ' 

Here the ion contribution (first term) is assumed as dispersionless, while we are 
interested in small wave vectors in comparison with the Brillouin-zone size. This 
term has the pole and the zero at the frequencies of the transverse and longitudinal 
phonons correspondingly. The natural phonon width iT na,t /2 is added to u; the 
high-frequency ion permittivity is denoted 6^. 

The second term in Eq. ([!]) is the electron contribution into permittivity for 
k <§C Pf, written |J with the help of the Boltzmann equation in approximation of 
the collision rate 7; the notation 

A p (k) = lo — kv + ry 

is introduced. The angular brackets 

, , J_ f, . 2dS f 



denote averaging over the Fermi surface with the density of states vq. For the isotropic 
case, the density of states on the Fermi surface u Q = m*pF/ir 2 and m* is the effective 
electron mass. 

Let us rewrite the electron contribution 



W L 1 - i( 7 /A p (fc)> 

where the Thomas- Fermi parameter k$ = 47ve 2 i'()/e 00 . 

For the isotropic Fermi surface we can carry out the integration 

(l/Ap(A;)> = — — ln^^ — . 

Ikvp 10 + — kvF 

4 



Separating the imaginary and real parts we have 



(1/A P (A;)) 



2kv, 



1 (uj + kv F ) 2 + 7 2 uj — kv F 
- In — — - + i arctan 

2 (uj — kvp) +7 7 



i arctan ■ 



uj + fcf ^ 
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(3) 



The imaginary part known as the Landau damping is pronounced at kv F > \uj + 



27 1 



e e (k, uj) = e OQ (ko/k) 2 (l + inuj /2kvi 
In the range 7 < uj — kv F « kv F , we have 



kv F » \uj + ry| 



k'n 



e e (k,uj) = Soo-g^ I 1 - 



UJ 



2kv F 



iln 



Ak 2 v 2 F 



— z- 
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(a; — kv F ) 2 + 7 2 a; — kv F 
For the small fcf^ << \uj + the expansion in kv F /(uj + ij) gives 



e e (k,uj) = e c 



1 



UJ 



pe 



3 / fcf F 
1 + 7 

5 V uj + 27 



(4) 



(5) 



(6) 



c<j(u; + 27) 

Here the ^-independent term represents the Drude conductivity. The electron plasma 
frequency is given by the integral over the Fermi surface 



UJ 



pe 



vdSf 



The limiting expressions are also valid for the arbitrary Fermi surface, but 

the constant v F has different values. In Eq. (Q), this is an average velocity on the 
belt v _L k, the velocity in the limiting point v || k in Eq. (^) , and the squared 
velocity averaged over the whole Fermi surface in Eq. (^). We do not pay attention 
the cases when the Fermi surface have flat or cylindrical pieces. 



III. FREQUENCY AND DAMPING OF PHONON PLASMON MODES 

The frequencies of the longitudinal phonon-plasmon modes are determined by 
the equation e(k,uj) = 0. In the absence of the electron and phonon collisions 
(7 = r nat = 0), we obtain with the help of Eqs. (|l|) and (|j) the biquadratic equation. 
It gives the frequencies of the coupled phonon-plasmon modes at k = 0: 



1 



1/2 
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These frequencies (related to the ujto) are shown in Fig. 1 (left panel) as functions 
of the electron concentration, namely, uj pe /uji . The upper dashed line begins at 
culo and tends to the electron plasma frequency uj pe . The lower frequency (solid line) 
starts as uj pe ujTo/^ho and then approaches c^to- In other words, observing in the 
optic range the longitudinal phonon mode and adding electrons, we see a transition of 
the longitudinal phonon frequency from uj lo to cjto This is a result of the Coulomb 
screening. 

Since both the collision rates 7 and T nat in the optic range are small in comparison 
with ujo, the damping of the phono n-plasmon modes can be added to their frequencies 
uj = uj± — iT±/2. Using Eqs. (H) and (|7|) we find at k = 0: 



r_ = [ 7 (a£ - ^o) + r nat (^ - <)]/(^ 2 - wi). 

The behavior of damping as a function of the electron concentration is shown in Fig. 
1 (right panel) for the ratio 7/r nat = 3. We see that the concentration dependence 
of widths corresponds with the type of modes: for the low electron concentration, 
the solid line represents mainly the electron plasma mode and the dashed line is 
associated with the longitudinal phonon. The character of modes reverses at the 
higher electron concentrations. 

Now let us consider the dispersion of the coupled phonon-plasmon modes. When 
the wave vector k increases, the upper mode (see Fig. 2 and Eq. @) approaches 
the asymptote uj = kvp\ 

, k\evl-ul Q ) 



u + = kvF < 1 + 2 exp 



k$(k 2 v 2 
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TO) 



— IJ. 



The lower mode drops into the domain kvp > u (dashed curve), where the Landau 
damping appears: 

2 _ k 2 ul Q + fcgu4 . nukkluli 



k 2 + k 2 Q 2v F (k 2 + kl) 2 ' v ' 

Here the ion plasma frequency ui pi = uj^q — ui\ - In the imaginary part, we should 
substitute instead of uj the frequency cu_ defined by the real part. Let us emphasize, 
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that the large dispersion arises under the conditions k VF 3> u and k < Pf- Then, at 
k ~ k Q , the frequency changes from uj^q to uj^o- Again this is the effect of screening. 

In the general case, the frequency and the damping of phonon-plasmon modes 
can be found in the numerical solution of the equation e(k, to) = 0. But of particular 
interest is the function — Im l/e(k, u), since it gives the intensity of the inelastic Ra- 
man or X-ray scattering, where ui and k have the sense of frequency and momentum 
transfers, respectively. The plots of the intensity, obtained with the help of Eqs. ([!]) 
- @ are shown in Figs. 3 - 5 for the cases of the large electron density (u pe > uto)- 
The peak at u/u^o = 2.6 -v- 2.7 corresponds to the plasmon excitation; notice, that 
u pe = k vp/\/3 for the quadratic electron dispersion. 

In the domain oj/u;to < 1 ; there is a phonon peak (see left panel of Fig. 1 and 
Eq. (0)). Because of the screening, its frequency is smaller than c^lo- F° r the small 
k (Fig. 3), the intensity of the phonon peak is 50 times less than the plasmon-peak 
intensity. For the larger k (see Fig. 4), the broad continuum u < kvp appears 
on the low-frequency side of the phonon peak. Here the dielectric constant e(k, u) 
has a noticeable imaginary part (|]) arising from the electron-hole excitations. The 
intensity of the phonon peak decreases, its line shape becomes asymmetric similar 
to the Fano resonance. Finally, the phonon peak is spread, being immersed in the 
electron continuum (see Fig. 5 and Eq. (H)). 

For the relative low electron density, u pe < cjtoj the intensities are shown in Figs. 
6-8. Now the phonon peak is approximately at uj/uo-^q = 1.5. That yields u = culoj 
since we set c^lo /^to = V%- There are also the plasmon peak and the electron- hole 
continuum. In Fig. 6 the plasmon peak is broader than the phonon peak, since we 
put 7/T nat = 3. The plasmon peak becomes more broader in Fig. 7 because of the 
neighboring electron continuum at u < kvp. It disappears completely in Fig. 8. 
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IV. CONCLUSION 



We have shown that the width of the longitudinal phonon-plasmon modes in- 
creases with increasing of the free-carrier concentration because their frequency ap- 
proaches the region kvp > ou, where the electron-hole excitations exist. The detailed 
experimental studies could help to clarify the comparative role of the different mech- 
anisms of the phonon shift and damping in the metallic and doped semiconducting 
materials under the influence of the electron concentration. 

The work was partially supported by the RFBR (project 01-02-16211). 
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Figure captions: 

Fig. 1. Frequencies (in units of ujto, left panel) and widths (in units of r nat , 
right panel) of the phonon-plasmon modes at k = in dependence of the free 
carrier concentration, namely, of the electron plasma frequency (in units of uito)- 
We set the ratio of the LO and TO frequencies in the absence of the free carriers 
^lo/Wo = \/2 and the ratio of the electron and phonon damping 7/T nat = 3. 
Fig. 2. The dispersion of phonon-plasmon modes for the metallic {uj pe > u;to, left 
panel) and semiconducting {u pe < o>ro, right panel) carrier concentration. The 
dashed straight lines separate the domain kv F > u, where the Landau damping 
exists; the dashed curves represent the strongly damped modes. 
Fig. 3-8. The inelastic light scattering intensity — Im l/e(k,u) as a function 
of the frequency transfer u for the momentum transfer k and the Thomas-Fermi 
parameter k (in the units of c^to/^f) shown in the figures. We set c^lo/Wo — V% 

r nat /cj TO = i(r 2 , 7/r nat = 3. 
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FIG. 2: 
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FIG. 6: 
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